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Abstract. It is known that quantum computers yield a speed-up for certain dis- 
crete problems. Here we want to know whether quantum computers are useful for 
continuous problems. We study the computation of the integral of functions from 
the classical Holder classes Fy a on [0, l] d and define 7 by 7 = (fc + a)/d. The 
known optimal orders for the complexity of deterministic and (general) randomized 
methods are 

comp(F^ ,Q , e) x e -1 / 7 

and 

comp random (F^V) x £-2/(1+27). 
For a quantum computer we prove 

comp^rM' V) x ^/(m) 

and 

COmpl^nt ( F &,a )£ ) < Ce -l/(l+7) ( log £ -l) 1/(1+7) . 

For restricted Monte Carlo (only coin tossing instead of general random numbers) 
we prove 

comp coin (F d fc 'V) < Ce- 2 ^ 1+2 ^ (loge- 1 ) 1 /^). 

To summarize the results one can say that 

• there is an exponential speed-up of quantum algorithms over deterministic 
(classical) algorithms, if 7 is small; 

• there is a (roughly) quadratic speed-up of quantum algorithms over randomized 
classical methods, if 7 is small. 



1. Introduction and Results 

1.1. Computation of the mean. Consider the following problem: Compute the 
mean 



1 n 



n 

i=l 



of n numbers Xj, where \xi\ < 1, up to some error < e < 1/2. The complexity of 
this problem depends on n and e and in the real number model we obtain 

(1) comp(ra, e) « n ■ (1 — e). 
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Here we consider the worst case setting, with the worst case cost and the worst 
case error. With randomized methods we can do much better, at least if n is large 
compared to e~ 2 . The cost is of the order 

(2) comp random (n, e) x min(n, e' 2 ). 

Now the error of a method is a random variable and the requirement is that its 
expectation is bounded by e. The statements ([!]) and (||) follow easily from well 
known upper and lower bounds. See, for example, Novak (1988). 

If we allow only random bits (restricted Monte Carlo methods, coin tossing) instead 
of arbitrary randomized methods then one gets the upper bound 

(3) comp com (n,e) < C ■ min(n, e~ 2 \ogn) 

which follows easily from (fj). 

For the results (1-3), and for all classical algorithms, we use the real number model 
of computation, with unit cost for each arithmetic operation. In addition we allow for 
(2) the instruction "choose a random number x G [0,1] "and for (3) the instruction 
"flip a coin" or "choose randomly {0, 1}" and also the cost of these instructions is 
one. See Novak (1995) for more details. 

A further improvement is possible by a quantum computer. The upper bound e~ l 
for the query complexity, defined by the number of times the real valued oracle is 
accessed to solve the problem, is proved in Brassard, H0yer, Mosca, Tapp (2000). See 
Theorem 12 of this paper. Using also the lower bounds from Nayak, Wu (1998) one 
can see that the exact order of this query complexity is 

(4) comp^(n, e) x min(n, e" 1 ). 

The bound (|j) is very important for the present paper. Grover (1998) states the 
upper bound e' 1 and says that it is optimal "up to polylogarithmic factors". 

If we consider, for the quantum computer, the bit number model then we need a 
slightly larger cost. H0yer (2000) proves the upper bound 

(5) comp quant (n,£) < C ■ mm(n , e~ l (logn + log loge -1 ), 

see Section || for details on the model of computation. The term n on the right side of 

(5) comes from the trivial classical algorithm. This term is n (logn + loge -1 ) in the 
classical bit number model. The output of a quantum algorithm is a random variable 
A(x,e), we always request that 

(6) \A(x,e) — S n (x)\ < e with probability at least 3/4. 

Of course we can run the algorithm several times and, taking the median from several 
measurements, we increase the probability of success. 
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1.2. Computation of integrals. How can we apply these results to the problem of 
numerical integration? Let us first consider the computation of the integral 



1(f) = / f(x)dx 

J[Q,l] d 

for functions from the Holder classes 

F2 = {/ : [0, l] d - R | H/IU < 1, \f(x) - f(y)\ < \\x - y\\ a }, 
< a < 1. Consider, for d = 1, the midpoint rule 



i=l x 



and for d > 1 the respective tensor product that uses n = £ d function values. By 
well known estimates for d = 1 together with the technique of Haber (1970, p. 489) 
we get the estimate 

(7) e(Q d n , F%) < C ■ d ■ n«l d 

for the worst case error of the product rule. To obtain e(Q d ll F£) ~ e, we have to take 

frd\ d/a 

(8) T»(i?, e )« — 



e J 

function values. We can now use the results from above to obtain upper bounds for 
the complexity of numerical integration. Using the (trivial) result (|l|) we get a bound 
for the (worst case) complexity of integration, 

f r d\ d/a 

(9) comp(F d V)< ^—j . 
With (g) we obtain 

(10) comp random (F d Q , e) < C ■ e~ 2 

and if we only allow random bits then we obtain, using (|3|) and (^|), 

(11) comp com (F d tt ,e) < C - e' 2 (\ogd + log^ 1 ). 

a 

In the same way we obtain for the quantum computer the upper bounds 

(12) comp^^^)^^^ 1 
and 

(13) comp quant (F d V) < C-e- l (\ogd + \oge- 1 ). 

a 
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Observe that all these bounds (|9|)-(|T3"D are just upper bounds which we get by a 
particular proof technique. Actually it is known that the order in (^|) is optimal, 

(14) comp(F«,e)^C dta e- d / a , 

while the upper bounds for Monte Carlo methods are not optimal, we have 

(15) comp random (F d a , e) « C d , a £ - 2d /^+ d ) 
and 

(16) comp coin (F d a , e) < C d , a £ - 2d/{2a+d) loge' 1 . 

For the proof of (|U|) see Heinrich (1993), Novak (1988), or Traub, Wasilkowski, 
Wozniakowski (1988). It is not difficult to show that ( pp follows from (0). Actually 
we will improve the exponent in the log-term slightly and prove such an upper bound 
with the factor (loge _1 ) 1 /( 1 + 2 «/4) > see (|31]). 



1.3. The problem and the results. Can the upper bounds ( [[4 ) and ( |13|) be im- 
proved, similarly as the upper bounds (|T0D and ([□])? What is the optimal rate of 
convergence (or the rate of the complexity) for numerical integration with a quantum 
computer? 

In this paper we answer this question for classes such as the F£ . We consider the 
more general Holder classes 

F k d > a = {} : [0, l] d ^ R I ||/ |U < 1, / G C\ |£>7(o0 - D*f(y)\ < \\x - y\\ a , VD 4 }, 



where D l runs through the set of all partial derivatives of order k and k G No, 
< a < 1. For k = we obtain F®' a = F£- It is convenient to use the notation 

k + a 

because this number is a good measure for the smoothness and appears in all the 
estimates. First of all, the optimal orders for deterministic and (general) randomized 
methods are known, see, e.g., Novak (1988). We have 

(17) comp^'V) -e- xh 



and 



comp random (i^V) ^ £~ 2/(1+27) . 



Therefore we have to study only the quantities comp quant , comp^, and comp com . 
For the upper bounds we use a technique called "variance reduction" in the literature 
on Monte Carlo methods. For the lower bound we use a decomposition technique of 
Bakhvalov, together with the lower bound of Nayak and Wu, see (|J). We obtain the 
following optimal rates of convergence. 
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Theorem 1. Define 7 = (k + a)/d, as above. 

(19) comp^(^^)^^ 1/(1+7) , 

(20) comp^i^'V) < C^e-WW (toge" 1 ) 1 ^), 

(21) comp coin (F*>V) < C d>k , a e' 2 '^ (loge" 1 ) 1 /^. 
To summarize the results one can say that 

• there is an exponential speed-up of quantum algorithms over deterministic (clas- 
sical) algorithms, if 7 is small; the multiplicative speed-up is roughly (1/e) 1 / 7 ; 

• there is a (roughly) quadratic speed-up of quantum algorithms over randomized 
classical methods, if 7 is small. 

1.4. Some comments. So far, most papers on quantum computing deal with dis- 
crete problems, such as factoring numbers or searching a database. Quantum comput- 
ing also helps for the continuous problems of numerical analysis or information-based 
complexity. Grover (1998) studies, among other things, the computation of the mean 
of finitely many real numbers. Related problems and algorithms are investigated also 
in the papers Boyer, Brassard, H0yer, Tapp (1998), Brassard, H0yer, Tapp (1998), 
Grover (1996), and Mosca (1999). The paper Nayak, Wu (1998) contains new lower 
bounds, while the recent paper Brassard, H0yer, Mosca, Tapp (2000) contains new 
upper bounds. 

Excellent surveys on quantum computing are Shor (1998) and Cleve, Ekert, Hen- 
derson, Macchiavello, Mosca (1999). Also the paper Abrams, Williams (1999) dis- 
cusses the computation of sums and integrals, "as long as the function is not patho- 
logical" . 

In numerical analysis and information-based complexity we usually assume the 
real number model with an oracle that gives function values, see Traub, Wasilkowski, 
Wozniakowski (1988) and, more formally, Novak (1995). For an expository account 
of continuous complexity on a classical computer see Traub, Werschulz (1998). Con- 
cerning the "allowed randomness" of the algorithms we may distinguish between 
three different cases. If a random number generator is available that can produce 
random u G [0, 1] according to the Lebesgue measure then we obtain the well known 
result (ITBT). Also the other extreme case, where no randomness is available, is well 



studied and we obtain the result (17). Hence we only have to consider the case of 
restricted Monte Carlo methods where coin tossing is allowed (and has unit cost), but 
not general random number generators. This case somehow corresponds to quantum 
computation where such a randomness can be easily realized. 

In Section |2| we give a little tutorial on quantum computation and present the 
search algorithm of Grover. We explain the model of computation and the cost of a 
quantum computation. Our proofs are contained in Section |3|. We add a section where 
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we discuss the computation of arbitrary bounded random variables by a quantum 
computer with a random number generator. Here we use a rather unrealistic model 
of computation because we assume that there is a random number generator without 
cost. 

2. The Model of Computation and the Search Algorithm of Grover 

In this section we describe the model of computation and the search algorithm of 
Grover (1996), see also Boyer et al. (1998). This section does not contain new results. 

Let Hi be a 2-dimensional Hilbert space over C and let eo and e± be two orthonor- 
mal vectors in Hi. The space H 1 represents a quantum bit, in the Dirac notation we 
have 

e = |0) and e x = |1) . 
For m G N quantum bits we use the 2 m -dimensional tensor product space 

H m = Hi®---®Hi 
with m factors. An orthonormal basis is given by the 2 m vectors 

be = e h <g> • • • <g> e im , 

where ij G {0, 1} and 

m 

(22) £ = Y,ij2 m ~ j , £ = 0,...,2 m -l. 

There are 2 m different b t and this corresponds to the 2 m different possibilities of an 
information that is given by m classical bits. The Dirac notation for bi is just \£), 
instead of e i± ® one finds |«! , z 2 ) or also \i<i). The formally different objects 
(ii, . . . ,i m ) and I or bi are often identified and called "classical state". 

One more piece of Dirac-notation is often used: \x) (y\ is a mapping, defined by 

(y\)(\z)):= (y,z)-\x). 

Here we write (y, z) — (y | z) for the scalar product. Therefore the projection P x on 
a normed vector x is written as \x) (x\. It is defined by y h- > (x, y) x. 
The Fourier series of x G iJ m is given by 

2 m -l 

(23) x = a (h,-,i m ) e h ® • • • ® e im = Y (3 e b e . 

ij£{o,i} e=o 

We are only interested in normed vectors, ||x|| = 1. All such vectors are called (pure) 
"quantum states" . For each quantum state there is a probability distribution on the 
classical states: the probability of £ is \j3e\ 2 . 
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A quantum algorithm starts with a classical state k G {0, . . . , 2 m — 1} which we 
identify with G H m . Then a number of unitary transformations Ui,...,U r are 
applied, the result is the quantum state 

x k = U r ...U l {b k ) 

and can be written in the form (p3]). Allowed are only those unitary transformations 
that are "efficient" in the sense that at most two quantum bits are changed. This 
means that, for example, Ui changes the first two bits and is of the form 

Ui{vx ® • • • ® v m ) = Ui(vi <g> v 2 ) <%> v 3 ■ ■ ■ ® v m , 

for some unitary Ui : C 4 — > C 4 . In the quantum bit number model, which we use 
for the numbers comp quant , one such unitary operation has cost one. The output of 
the algorithm, given by a final measurement, is a classical state I G {0, . . . , 2 m — 1}, 
or certain bits of i. The probability of i is |/^| 2 , where is the respective Fourier 
coefficient of x k = U r . . . Uiipk). We say that a quantum algorithm computes a given 
function if the probability of a correct output is at least 3/4. 

Now we describe the search problem and quantum computations with an oracle. 
Let m G N and X m = {0, l,...,2 m - 1}. Assume that / : X m -> {0,1} is an 
arbitrary mapping which, of course, can be identified with a subset of X m . We define 
a corresponding unitary mapping Sf on H m by 

S f (b e ) = -b t if = ! and S f {b t ) = b & if = 0. 
We also put So = S/ for /(£) = <5 0> £. A 6/acA; box or oracle Qf for / is defined on 
# m +i by 

Qf{bt® ej) = &^ ® e i+/( £). 
Here the plus sign in e i+ /^) means addition modulo 2, also called exclusive or. Then 
one can easily show that 

S f (bi) ® e = Q f PQf(b e <g> e ), 

and therefore the oracle can be used to compute function values of Sf. Here P is 
defined by 

P(b t ®ei) = (-lybiQet. 
One can even simulate Sf with Qf using only one application of Qf\ simply apply 
Qf on (b e ® (e - ei)). 

^4 search problem is defined as follows. Let F m be the set of all fi : X m — > {0, 1} 
with fe(j) — 1 iff j = £. Of course we may identify the sets F m and X m , and to each 
£ G A m or G F m there is exactly one bg. The problem is to find I if / = ft is given 
by the oracle Q/. 

The algorithm of Grover works as follows. First we define the Walsh-Hadamard 
transform W\ : H\—> H\ by 

W 1 (ei) = ^(eo + (-l)*e 1 ) 
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and W m = W\ ® ■ ■ ■ <S> W±. Now the algorithm is defined by 

(24) {-W m S W m S f ) k (W m (b )), 

where k has the order 2 m//2 . It is shown in Boyer, Brassard, H0yer, Tapp (1998) that 
k « 7r 2 m / 2-2 is a very good choice that leads to a high probability of success. 

The cost, in the quantum bit number model, of every iteration in fl24]) is about m, 
the total cost (to find the element I with high probability) is about m ■ 2 m l 2 or, with 
iV = 2 m , about y/N logiV. Formally the algorithm is slightly different because Sf is 
only given by the oracle Qf, hence we work with m + 1 instead of m quantum bits. 
We assume that an application of Qf has unit cost. 

For the problem "compute the mean of xi, . . . ,x n " we assume, when we consider 
the quantum bit number model and the numbers comp quant , that there is a Boolean 
oracle Q that gives the jth bit (digit) of Xj. For the quantum query complexity 
we allow real valued oracles and only count the number of oracle calls. In the case 
< if < 1 a real number quantum oracle has the form 

Q(b t (g> e ) =bi>® <Jx~t e + be® y/l — xe e\. 

For all classical algorithms we allow real- valued oracles for the X£. 

3. Proof of Theorem 1 

We fix a space F^' a . First we prove the upper bounds for quantum computers. We 
use an algorithm of the form 

(25) A(f)=I(P n f)+Q%(f-P n f). 

Here I(P n f) is the integral of a function P n f and Q d N is the d- dimensional midpoint 
rule, as in Section 1.2, which we apply to (/ — P n f). We now explain the operators P n 
and Qfj. By P n we mean a projection operator by interpolation, one can use piecewise 
polynomials, which uses n function values and gives an order 

\\f -Pnf \\oo~ n~\ 

for / G F^' a . It is well known that this is the optimal order of convergence, see Novak 
(1988). One evaluation of (/ — P n f) can be implemented at a constant cost (where 
the constant depends on d and k, but not on n). 

By Q d N we mean that we do not really apply the midpoint rule Q d N . Instead we 
evaluate this midpoint rule by a quantum computer up to some error e\ ■ rT 1 with 
the cost 

COSt^y^O^C-^ 1 

or 

cost quant (JV,ei) < C-^logiV + loglog^ 1 ), 
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respectively. The error of this method is bounded by 

(26) e(f)<C.N^+n-^e 1 , f e F%' a , 

which is the sum of the error by discretization (the integral being replaced by Q d N ) 
and the error made by the approximate evaluation of Q d N . This error bound is valid 
for all (5 < a/d, if k — 0, and /3 = 1/d, if k > 0. In addition we have to assume that 
N is at least of the order n. To simplify the presentation we use a (3 which is always 
smaller than 1. The complete cost of the method is bounded by 

(27) cost^(A^i)<C7> + 5r 1 ) 
or 

(28) cost quaBt (JV,^) < C-n + C-e^QogiV + loglog^ 1 ), 
respectively. In the query-complexity case we simply choose n « e^ 1 and 

AT' 3 « n" 7 • e a 

and observe that we may apply fl2~6|) , because of (3 < 1. We obtain a cost of the order 
n and an error of the order n~ 7_1 and so obtain the upper bound in (|T^). 
In the quantum bit number model we take 



again with 

Then we get the cost bound 
and the error bound 
We obtain 



n « e x 1 \oge l 1 , 
N' 13 « n~ 7 • ei. 

cer 1 log^r 1 

e(/)<Ce 7+1 (log^ 1 )" 7 - 
e (/)-V(7+D > C£ -i (bg^ 1 ) 7 /^ 1 ) 



and therefore 

comp quant (F d fc ' 0, ,e) < Ce" 1 /^ (loge" 1 ) 1 - 7 ^^. 

Now we prove the lower bound in (p!9|). The space i 7 ^'" contains n x e 1 1/7 functions 
/i, /2, . . . , f n with disjoint supports such that 

• /[o,i] d f* dx = £ i +1/7 and 

• EIli Ai/i 6 F d fc ' Q if |A,|<1, 
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see Novak (1988, p. 35). Consider now the following problem. Compute the mean 
value of the integrals L ^ d Aj/j dx, where |Aj| < 1, up to some error e 2 . We can apply 
the lower bound of Nayak, Wu (1999), see (^), to obtain the lower bound 

cost >Cmin(£~ 1/7 , e\ +lh e 2 l ). 

Of course we put e 2 = ej +2 ^ 7 and obtain 

cost > Ce1 lh . 

What we estimated was the cost to compute the mean value. Since Jj Q x , d Ym=i ^ifi dx 
is actually the sum, the error is to be multiplied by n, hence s = e 2 n ~ e\ +1 ^. Since 
the cost to obtain this error is at least of the order e^ 1 ^ we obtain 

oomp^i^e)^" 1 /^. 

We finally prove the upper bound for restricted Monte Carlo. Instead of the arbi- 
trary random numbers of a (general) Monte Carlo method we can only use random 
bits or coin tossing. We use a discretized version of a well known variance reduction 
technique and write the method in the form 

(29) A(f) = I(P n f) + Q d N (f-P n f), 

and only the meaning of Q d N is different from By Q d N we mean that we do not 



really apply the midpoint rule Q d N . Instead we evaluate this midpoint rule by the 
classical Monte Carlo method, again up to some error E\ ■ n -7 with the cost 

cost coin (iV, £i) < C • log N + n. 

The bound e± 2 is the classical Monte Carlo bound, the factor log N comes in because 
we need (about) logiV random bits to select one node from the possible iV nodes. 
The error of this method is bounded by 

e(f) < C ■ N-P + ■ e 1} f E F^' a , 

see (^6|). Again we use a (3 which is always smaller than 1. We put ps* rT 1 E\ and 
n ps e\ 2 logej -1 and obtain 

cost coin (iV,£i) < C -e\ 2 loger 1 

and 

e(f) < Cn^e 1 < C e} +27 (log^V 7 

and therefore 



QUANTUM COMPLEXITY OF INTEGRATION 11 

4. A Remark on Randomized Quantum Algorithms 

We present an algorithm to compute the expectation of arbitrary bounded random 
variables up to some error e > 0. A classical randomized method needs time e~ 2 . The 
proposed algorithm uses the algorithm of Grover and a random generator. We assume 
that the random generator is for free and this is certainly not a realistic assumption. 
Let (X,B,m) be a probability space and let 

F = {/ : X — > R | / is measurable and ||/||oo — !}• 

We want to compute the integral of a function (or the expectation of a random 
variable) 

Hf) = [ f(x)dm(x), 
Jx 

for / G F. With (||) one obtains the complexity bound 

(30) comp- n a d n r(F, e) < C e~ x \oge~\ 

with a constant C that does not depend on the particular space (X, B, m). If we use 
the quantum query model, together with a free random generator, then we get in the 
same way the upper bound C e^ 1 . We prove fl3*0|). 

To compute an approximation A(f,e) of /(/) for / G F and < e < 1/2 (with 
the understanding that @ should be true) we proceed as follows: 

• First we randomly select X\, . . . , x n G X using the random generator, where 

(31) n=\72e- 2 ]. 
If we putQ 

1 n 
n L — ' 

i=l 

then it follows from Chebyshev's inequality that 

(32) |/(/) - Q n (f)\ < - with probability at least 7/8. 

3 

• We assume that an oracle is available for the computation of f(xi), up to a 
maximal error of e/3. On input i the oracle gives f(xi) such that 

|/(x 4 )-/(*i)|<|. 



x Of course Q„ is simply the classical Monte Carlo algorithm for the computation of the integral. 
Here we only define the points x\, . . . , x n , we do not compute Q n (f). This expression is only used 
for the analysis and the intuition. 
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• Quantum computation. Let 

i n 

Qn(/) = -^fe). 
i=l 

Then we know \Q n {f) — Qn(f)\ < ^/3 and therefore 
~ 2 

l-^(/) — Qn(f)\ < q £ with probability at least 7/8. 
Hence we compute an approximation A(f,e) of Q n (f) such that 

(33) \A(f, e) - Q n (f)\ < - with probability at least 7/8 

3 

and get (^|). With fl5|) and (31) one obtains the bound fl3~0|) for the complexity 
of the problem. 
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